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Investigation of variable star candidates in the globular cluster 
NGC 5024 (M53) 
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ABSTRACT 

We have performed a ca reful investigation of the 74 candidate variable stars presented by 
ISafonova & Stalin ( 201lh. For thi s pu rpose we used our data base o f imaging and light curves 
from Arellano Ferro et al.l (1201 lb and lArellano Ferro et al .1 (120121) . We find that two candi 



dates a re known variable stars, eight cand idates were discovered fir st bv lArellano Ferro et al.l 
( 20 1 1 ) but would not have been known to S afonova & Stalinl (1201 ll) at the time of their paper 
submission, while four candidates are new variables. Three of the new variables are SX Phe 
type and one is a semi-regular red giant variable (SR type). We also tentatively confirm the 
presence of true variability in two other candidates and we are unable to investigate another 
four candidates because they are not in our data base. However, we find that the remaining 
54 candidate variable stars are spurious detections where systematic trends in the light curves 
have been mistaken for true variability. We believe that the erroneous detections are caused 
by the adoption of a very low detection threshold used to identify these candidates. 

Key words: stars: variables - globular clusters: inidividual: NGC 5024 



1 INTRODUCTION 

Performing a census of variable stars in a globular cluster is an im- 
portant task because variable stars, such as RR Lyrae and SX Phe, 
provide a wealth of information about the host cluster. Fourier de- 
composition of the RR Lyrae light curves can be used to estimate 
their metallicity and absolute magnitudes, from which estimates of 
the metallicity and distance to the cluster may be derived. A fre- 
quency analysis of the light curves of the SX Phe stars combined 
with the knowledge of their period-luminosity (P-L) relation may 
also be used to determine a cluster distance, and conversely, the 
results may be used to investigate how the SX P he P-L relation 
depends on metallicity ( Cohen & Sarajedi nil2012l) . 

When performing precise differential time-series photometry, 
essential for a variable star census, there is always the problem of 
systematic photometric errors that correlate with image and object 
properties such as point-spread function full-width half-maximum 
(PSF FWHM ), spatial coordinates, local sky background, etc. ("red 
noise " - see IPont. Zucker & Oueloj [2006: Bra mich & Freudlind 
l20l2h . Systematic trends may sometimes be mistaken for real vari- 
ability, since image and object properties may vary smoothly over 
time, mimicking variability on various timescales. Therefore it can 
be quite easy to contaminate the catalogue of true variable stars in a 
globular cluster with non- variable stars whose light curves exhibit 
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systematic trends, and consequently these stars may be accidentally 
used as part of the sample of variable stars to derive the cluster 
properties, introducing unwanted errors into these estimates. One 
way to avoid contamination is to set the detection threshold for 
variables to a level that has a very small probability of admitting 
false detections even when systematic trends are present. The de- 
tection threshold may be determined empirically from an analysis 
of the light curves of the general ensemble of stars which are usu- 
ally dominated by true non-variable stars. 

In a series of papers (e.g. lAre llano Ferro et al. 



lArellano Ferro. Giridhar & Bramichl l2010t iBramich et al.l 1201 ill 
we have been performing a variable star census for a range of 
Galactic globular clusters of different metallicities. We have tar- 
geted clusters that are generally under-studied with regards to time- 
series photometry in the CCD era, and we aim to detect all the vari- 
ables in the cluster field-of-view down to close to our magnitude 
limit. To avoid confusion over the variable star identifications for 
future researchers, we provide detailed finding charts and accurate 
astrometry in addition to the light curves of the variable stars that 
we find. 

Specific ally, we recently perfo rmed this task for 
NGC 5024 ^Arellano Ferro et all I20T1I - from now on All; 
Arella no Ferro et al .1120121 - from now on A12). We used the 2.0 m 
telescope of the Indian Astronomical Observatory (IAO), Hanle, 
India, equipped with an imaging camera with a pixel scale of 0.296 
arcsec pix~' and a field-of-view (FOV) of ~ 10. 1 x 10.1 arcmin 2 , to 
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obtain V and / time-series photometry with ~300 epochs in each 
waveband observed during 18 nights spread out over ~2 years. 
Our exposure times ranged from 25 s to 600 s in both filters, with 
the longer exposure times generally employed during the earlier 
observation nights. We reduced the image data using the technique 
of difference imag e analysis (DI A) and specifically we used the 
DanDIA software (Bramich 2008) for this purpose. We clarified 
the variable status, periods, and ephemerides of the previously 
known variables, including correcting the misidentification of 
three variables, and we presented the discovery of two new RR 
Lyrae stars and 13 new SX Phe stars. 

More recently. ISafonova & Stalinl d201lh (from now on SSI 1) 
used the same telescope and camera to obtain ~7.2 hours of con- 
tinuous ^-filter time-series photometry (101 epochs) from a single 
night and using a single exposure time of 100 s. We therefore ex- 
pect that the SS11 image properties are very similar to those of 
the All and A12 images (e.g. cluster positioning, seeing, number 
of saturated stars, etc.). SS11 also reduced their image data using 
DanDIA, although the parameters used for the reduction process 
were not necessarily the same as those used in Al 1 and A12. SS 1 1 
claim the detection of 74 new variable stars from the analysis of 
their data. The purpose of this paper is to investigate these new can- 
didate variable stars in order to confirm their variable status using 
our more extensive photometric data base. We also take this op- 
portunity to publish some corrections of mistakes in the literature 
regarding variable star identification and astrometric coordinates in 
Appendix A. 



2 INVESTIGATION OF THE SS11 VARIABLE STAR 
CANDIDATES IN NGC 5024 

A total of 74 stars are presented as candidate variables in SS 1 1 ; 
namely, one standard star exhibiting suspected eclipses, 14 RR 
Lyrae-type variables, 10 eclipsing binaries, 22 SX Phe-type vari- 
ables, and 27 other unclassified variables. In this Section, we inves- 
tigate each set of candidate variables in turn. 



2.1 Variability In Standard Stars 

SS 1 1 note that six Stetson standard stars dStetsonl200(j|) exhibit spu- 
rious variability with the same pattern in their light curves, which 
they explain as being caused by a deficiency in the data reduction 
process. However, star S240 is also suspected of showing the sig- 
nature of an Algol-type detached eclipsing binary with two sus- 
pected eclipses of depths ~0.05 and ~0.1 mag separated in time 
by ~2.5 h, although ultimately SS 1 1 classify the variability seen in 
their light curve as spurious (see their Table 5). 

We plot our V light curv^]for S240 in Figure Q]using the same 
scale on the magnitude axis as that used in SS11. Our light curve 
shows some systematic photometric trends at the level of ~0.01- 
0.02 mag and we achieve a root-mean-square (RMS) magnitude 
deviation in the light curve of ~9 mmag. Given that our light curve 
includes multiple sections of continuous time-series photometry of 
durations greater than 4 h, we should be able to detect an eclipsing 
signal with period ~2.5 h and amplitude ~10 times our light curve 
noise. However, we detect no such signal, and we conclude that the 



The / light curve does not exist in our data for this star because it lies on 
a bad column in the / filter reference image. 



eclipse-like signatures seen in the SSI 1 light curve of star S240 are 
artefacts of the data reduction process. 

2.2 RR Lyrae-Type Variable Candidates 

As described in iBramich et al] d201ll) . the DanDIA software mea- 
sures the differential flux of each star in the difference images 
by optimally scaling the appropriate point-spread function (PSF) 
model at the known star position (as determined from the refer- 
ence image) while simultaneously fitting for a local background. 
The software stores the PSF model as a circular image stamp of 
radius four times the FWHM and it follows that for star pairs in 
close spatial proximity their corresponding PSF models may over- 
lap. However, the DanDIA software only performs an independent 
fit to each star in a difference image and consequently the differen- 
tial flux from a real variable may systematically influence the re- 
sults of the PSF fit to any nearby stars, with the effect being largest 
for the closest and faintest nearby stars. An easy way to spot if the 
flux variations from a nearby variable star have influenced the pho- 
tometry of an otherwise constant star is to see if the light curves of 
both stars correlate, or if they have the same period and phase. 

SS11 present 14 candidate variables that are suspected of 
showing RR Lyrae-type light curve variations. However, we find 
that candidates RR1, RR2, RR8, RR12, and RR14 lie within ~7.6, 
3.8, 16.9, 11.4, and 20.2 pix, respectively, of the known RR Lyrae 
variables VI 1, V18, V40, VI, and V58, respectively. Given that the 
typical PSF FWHM of the SSI 1 observations is ~5 pix (~1.5 arc- 
sec), these candidate variables have PSFs that overlap substantially 
with the nearby RR Lyrae star PSFs. Furthermore, these candidate 
variables are fainter than the known RR Lyrae variables (by ~1.5- 
2.5 mag) and therefore we believe that the RR Lyrae-like variations 
that are seen in the SS11 light curves of RR1, RR2, RR8, RR12, 
and RR14 are due to the systematic influence on the light curves 
from the nearby RR Lyrae stars. 

Using the celestial coordinates from SS11 of the remaining 
RR Lyrae-type variable candidates, we have identified these stars 
in our V and / reference images and found the corresponding light 
curves in our data base. Candidates RR4 and RR5 are very faint 
in our reference images and we do not have the corresponding 
light curves. We cannot therefore comment on the variations seen 
in these stars by SS11. For RR3, RR6, RR7, RR9, RR10, RR11, 
and RR13, we plot our light curves in Figure Q] using black and 
red points for the V and / photometric measurements, respectively, 
and using the same scale on the magnitude axis for each candidate 
variable as that used in SSI 1 to enable a direct comparison. The / 
light curves have been offset in magnitude by an amount that al- 
lows a clear comparison of the light curves between the two filters. 
Note that the scatter in our light curves increases with time simply 
because our adopted exposure times were systematically reduced 
throughout our observing campaign. 

True photometric variability generally correlates between dif- 
ferent filters whereas systematic trends in the photometry do not 
necessarily correlate but they may do. Furthermore, systematic 
trends tend to exhibit similar temporal patterns between stars that 
are relatively close to each other in an image. Candidates RR1 1 and 
RR13 are within 50 pix of each other and their V light curves show 
similar trends with time, which is also the case for the R lightcurves 
of these stars displayed in Figure 11 of SS 1 1 . Also, these candidate 
variables are very close to a star that is highly saturated in both of 
our reference images which has a negative impact on the quality 
of the difference images in this region, explaining the particularly 
large amplitude of the trends in our light curves for these stars. 
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Figure 1. Light curves from All and A12 for selected candidate variables from SS 1 1 . V and / filter light curves are plotted in black and red, respectively, with 
mean photometric uncertainties per data point on each night plotted at the beginning of a night as vertical bars. The magnitude scale for each candidate has 
been matched to the scale used in S S 1 1, and the / filter light curves have been offset in magnitude by an arbitrary amount for clarity. 
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Later on in Section 1231 we also conclude that the variable candi- 
dates VC18, VC19, and VC20 exhibit spurious variability due to 
their proximity to the same highly saturated star. For these reasons 
we do not believe that SS11 have detected true variability in stars 
RR11 and RR13. 

Our light curves for candidate variables RR7 and RR10 rule 
out the existence of the ~0.2-0.3 mag amplitude variations seen in 
the SS11 light curves of these stars. Inspection of our light curves 
for RR3, RR6, and RR9 reveals that the trends present in our light 
curves do not generally correlate between filters (e.g. see the night 
of 17/04/2009), and certainly the trends do not reach the ~0.5 and 
0.4 mag amplitudes seen in the SS 1 1 light curves for RR6 and RR9, 
respectively. We therefore believe that the larger-amplitude varia- 
tions seen by SSI 1 and the smaller-amplitude variations seen in our 
light curves for these five candidates are due to trends introduced 
during the reduction process. 

Finally, in this section we comment that all of the RR Lyrae- 
type variable candidates RR1-RR14 as reported by SSI 1 are fainter 
(R ~ 18 -20 mag) than the horizontal branch of NGC 5024 (R ~ 
17.2 mag), which implies that if they are RR Lyrae stars, then they 
must lie behind the cluster relative to our Sun, and therefore they 
cannot be cluster members. It is highly unlikely that 14 Galac- 
tic RR Lyrae stars not associated with NGC 5024 lie within this 
~10x 10 arcmin 2 field of view. This con clusion is supported by the 
fact that over an area of ~ 10000 deg 2 . ISesaretal.N201ll) found 
on average only ~0.5 RR Lyrae stars per square degree down to 
R ~ 18 mag. If any of these stars are eventually found to be true 
variables, then they are most likely to be of the SX Phe type. 

2.3 Eclipsing Binary Variable Candidates 

SS11 present 10 candidate variables that they classify as short- 
period contact eclipsing binaries. They determine a period for each 
candidate and present the phase-folded light curves. Using the ce- 
lestial coordinates from SS11, we have found the corresponding V 
and / light curves in our data base. 

Candidates Wl and W8 lie within ~14.2 and 15.4 pix of the 
known variables V33 (RR Lyrae type) and V68 (SR type - red giant; 
highly saturated in our data), and for the same reasons as described 
in Section [2~2l we do not believe that SSI 1 have detected true vari- 
ability in these stars. However, candidate W9 is a true variable star, 
but it is not an eclipsing binary. Instead it is the SX Phe type vari- 
able V104 first detected by All. SS11 derive a tentative period of 
~0.15 d for V104 consistent with the much more precise period of 
~0. 14756 d from All. 

The periods determined by SSI 1 for any of their variable can- 
didates are not appropriate for folding our light curves over the 
full time span of our observations (~2 years) since they are de- 
rived from a time base of only ~7.2 hours. Therefore, whenever we 
choose to plot phase-folded versions of our light curves using the 
periods from SS 1 1, we only plot the photometric data from the first 
three nights of observation (17-19th April 2009) since these nights 
are consecutive and the data have the best photometric precision 
in our time-series because they correspond to the longest exposure 
times of 600 s in V and 420 s in /. 

In Figure [2] we plot our light curves for candidates W2-W6, 
Wll, and W13, phase-folded using the periods from SS11. We do 
not have an / light curve for candidate W13 because this star lies 
on a bad column in the / filter reference image. SSI 1 list the am- 
plitudes of the light curve variations in the R band as ~1.22, 0.49, 
0.75, 0.30, 0.60, 0.54, and 1.06 mag for the candidates W2, W3, 
W4, W5, W6, Wll, and W13, respectively, and these variations, 



if true, will have similar amplitudes in the V and / bands. Inspec- 
tion of our plotted light curves reveals that we can fully rule out the 
presence of these variations in the candidates W2, W4, W5, W6, 
and W13. For W3, our light curves show a relatively large scatter, 
but it is clear that there are no coherent variations at the SS 1 1 pe- 
riod. We therefore believe that the variations seen by SSI 1 are due 
to trends introduced during the reduction process. 

For the variable candidate Wll, our phase-folded light curve 
in Figure [2] shows a hint of coherent variations at the SS11 pe- 
riod that match the shape of the SS11 phase-folded light curve, 
and therefore we tentatively confirm the variable nature of this star. 
However, the scatter in our light curve is sufficiently large that we 
believe that better precision follow-up observations are required 
before this candidate variable is assigned a proper V-number. We 
note that Wl 1 has V sa 20.306 mag and V - 1 0.763 mag, which 
places it outside the blue straggler region as marked in the colour- 
magnitude diagram (CMD; Figure 4 in All). We also find that its 
period and mean magnitude are inconsistent with the P-L relation 
for SX Phe cluster members (see later in Figure|6], since the star is 
much too faint for such a long perioc0. We therefore cannot specu- 
late on the type of variability exhibited by Wl 1. 

To quantify our conclusions about the candidate variables 
plotted in Figure [2] we calculate the improvement in chi-squared 
A% 2 when fitting a sine-curve compared to a constant magnitude. 
Under the null hypothesis that a light curve is not variable, the Ax 2 
statistic follows a chi-squared distribution with two degrees of free- 
dom. We set our threshold for rejection of the null hypothesis at 1%, 
which is equivalent to A^ 2 > 9.21. Candidate Wll has A^ 2 « 9.33 
and 9.39 for the V and / light curves, respectively, supporting our 
tentative conclusion that it is variable. We plot the fitted constant 
and sine-curve models in the corresponding panel of Figure [2] as 
the dashed and continuous curves, respectively. For the remaining 
candidates shown in Figure[2] the maximum value of A^ 2 is ~4.42, 
meaning that we do not detect any sinusoidal variability in our data 
at the SS11 periods. 

2.4 SX Phe-Type Variable Candidates 

SS11 present 22 candidate variables that they classify as "poten- 
tial SX Phe stars". They determine a period for each candidate and 
present the phase-folded light curves. Again, using the celestial co- 
ordinates from SS11, we have found the corresponding V and / 
light curves in our data base, except for the candidate SX13 where 
we cannot find a star at the SS 1 1 coordinates in our V and / refer- 
ence images. 

Candidate SX6 is the known SX P he variable star V78 first 
discovered by iDekanv & Kovac j 12009]) . Candidates SX2, SX3, 
SX9.1, SX11, SX14, SX19, and SX22 are the SX Phe variable 
stars V96, V97, V99, V101, V102, V103, and V105 first discov- 
ered by All. Candidate SX9, as noted by SS11, is only ~4.2 pix 
from the true variable V99, and its SS 1 1 photometry has obviously 
been influenced by the flux variations from V99 because SS11 de- 
rive very similar periods and the same epoch for both stars. Also, 
as acknowledged by SS11, candidate SX25 is the known S X Phe 
variable star V80 first discovered bv lDekanv & Kovacsl d2009l) . 

In Figure [3] we plot our light curves for the candidates SX4, 
SX7, SX8, SX12, SX15-17, SX20, SX21, SX23, and SX24, phase- 
folded using the periods from SS 1 1 . We do not have an / light curve 
for the candidate SX15 due to a failure of the DanDIA software to 

2 Wl 1 is not plotted in Figure [6] 
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Figure 2. Light curves of the first three nights of observations from Al 1 and A12 for selected candidate variables from SS11. The light curves have been 
phase-folded using the periods derived by SSI 1. V and / filter light curves are plotted in the left and right-hand panels, respectively. The magnitude scale for 
each candidate has been matched to the scale used in SSI 1. For the candidate Wl 1 showing true variability, fitted constant and sine-curve models are plotted 
as dashed and continuous curves, respectively. 



measure the reference flux of this star in the / filter reference image. 
Applying the same chi-squared analysis to these folded light curves 
as that described in Section|23] we find that SX4 has A# 2 14.43 
and 10.22, SX12 has Ax 2 » 8. 18 and 10.86, SX16 has A^ 2 « 23.05 
and 3.79, and SX21 has A^ 2 « 34.67 and 1 1.59, all for the V and 
/ filters, respectively. We therefore confirm the variability at the 
SS11 periods for the candidates SX4, SX12, SX16, and SX21, al- 
though we consider the candidate SX12 as border-line because the 
two A^ 2 values for the light curve lie only just above and below 
the detection threshold. We plot the fitted constant and sine-curve 
models in the corresponding panels of Figure [3]as the dashed and 
continuous curves, respectively. 



For the remaining candidates SX7, SX8, SX15, SX17, SX20, 
SX23, and SX24, the maximum value of A^ 2 is ~3.51, and we 
conclude that our data do not show any variability at the SSI 1 pe- 
riods. Also, if the variations detected by SSI 1 in their light curves 
for these stars are due to true variability, then we should have been 
able to detect them since their measured amplitudes are consider- 
ably larger than the scatter in our light curves in all cases. Therefore 
we believe that the variations seen in the SS 1 1 light curves of these 
candidates are due to trends introduced during the reduction pro- 
cess. 



2.5 Unclassified Variable Candidates 

SS11 present 27 candidate variables "that show definite light vari- 
ations" but which they are unable to classify. Using the celestial 
coordinates from SS11, we have found the corresponding V and / 
light curves in our data base, except for candidate VC9 where the 
DanDIA pipeline has failed to measure the reference flux of this star 
in both the V and / filter reference images. 

Candidates VC1, VC2, VC7, VC11, VC15, VC16, VC17, 
VC21, VC22, and VC24 lie within ~3.7, 4.0, 3.7, 20.4, 20.2, 3.5, 
13.3, 4.4, 4.7, and 4.3 pix, respectively, of the known RR Lyrae 
variables V33, V42, V64, V46, V72, V53, V63, V45, V9, and V51, 
respectively. Inspection of the SS11 light curves for these candi- 
dates reveals that they mimic segments of RR Lyrae light curves, 
and we believe that the variability detected by SS 1 1 for these stars 
is due to the systematic influence on the light curves from the 
nearby RR Lyrae stars. 

Candidates VC18, VC19, and VC20 all lie within ~10 pix of 
the same highly saturated star that is also in the vicinity of the pre- 
viously discussed candidates RR1 1 and RR13 (see Section [2~2t . We 
believe that the poor quality of the difference images in such close 
proximity to a highly saturated star is the cause of the variations 
observed in the SS 1 1 light curves for these stars. 

In Figure [4] we plot our light curves for VC3, VC5, VC6, 
VC8, VC10, VC12-14, and VC25-28 using the same scale on the 
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Figure 3. Light curves of the first three nights of observations from Al 1 and A12 for selected candidate variables from SS11. The light curves have been 
phase-folded using the periods derived by SS 1 1 . V and / filter light curves are plotted in the left and right-hand panels, respectively. The magnitude scale for 
each candidate has been matched to the scale used in SSI 1. For the candidates showing true variability, fitted constant and sine-curve models are plotted as 
dashed and continuous curves, respectively. 



magnitude axis as that used in SS11, and in Figure Q] we plot out- 
light curve for V8CL7, which was flagged as a suspected variable 
by SS11 during their investigation into the identification of V80. 
As expected, all of our light curves show systematic trends at some 
level (usually the few per cent level). However, for all of the can- 
didates except VC27 and VC28, the amplitudes of the variations in 
our light curves are considerably smaller than the amplitudes of the 
variations seen in the SS 1 1 light curves. Furthermore, the variations 



seen in our light curves do not generally correlate between wave- 
bands on a reasonable fraction of nights. We also note that VC25 
and VC26 are close to a saturated star. Hence, on the basis of our 
data, we do not confirm the variability of the stars VC3, VC5, VC6, 
VC8, VC10, VC12-14, VC25, VC26, and V80.7, and we believe 
that the variations seen by SS11 are most likely the signature of 
systematic errors introduced during the reduction process . 

For candidate VC27, SS11 present a light curve with varia- 
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Figure 4. Light curves from Al 1 and A12 for selected candidate variables from SSI 1. V and / filter light curves are plotted in black and red, respectively, with 
mean photometric uncertainties per data point on each night plotted at the beginning of a night as vertical bars. The magnitude scale for each candidate has 
been matched to the scale used in SS 1 1 , and the / filter light curves have been offset in magnitude by an arbitrary amount for clarity. 
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V106 



V107 





V108 



V109 



Figure 5. Finding charts constructed from our V reference image; north is 
up and east is to the right. Each image stamp is of size 23.7x23.7 arcsec 2 . 
Each confirmed variable lies at the centre of its corresponding image stamp 
and is marked by a cross-hair. 



tions of amplitude ~0.15 mag, and our light curve is not of suffi- 
cient precision to definitively rule out such variations. However, the 
variations that are seen in our light curve for VC27 do not always 
correlate between wavebands (e.g. see the night of 19/04/2009). For 
candidate VC28, our light curve on any single night does not show 
the ~0.3 mag variations seen in the SSI 1 light curve, and again the 
variations do not always correlate between wavebands (e.g. see the 
night of 17/04/2009). Hence, we also do not confirm the variability 
of the stars VC27 and VC28. 



3 NEW VARIABLES IN NGC 5024 



In Section [2T4| we confirm the variable nature of the candidates 
SX4, SX16, and SX21 from SS11, and consequently, we assign 
them the new variable numbers V106, V107, and V108, respec- 
tively. For future ease of identification, we provide the finding 
charts for these stars in Figure [5] We perform a frequency anal- 
ysis of the light curv es of these new variables using PERI0D04 
(Lenz & Bregedl2005l) , the results of which we report in Table Q] 
along with their mean magnitudes, colours, and celestial coordi- 
nates. 

We only detect one dominant frequency for each variable, 
which we use to plot the stars in Figure[6]as solid red circles. This 
figure displays the P-L relation for the SX Phe stars in NGC 5024 
and it is an updated version of Figure 12 from Al 1. VI 06 and V 108 
follow the fundamental mode P-L relation for the cluster and they 
also lie in the blue straggler region of the CMD (Figure 4 of Al 1). 
Hence we classify V106 and V108 as SX Phe stars pulsating in 
the fundamental (F) mode. V107 lies above the P-L relations for 
the fundamental, first, and second overtone pulsators in the cluster, 
and it also lies close to but outside of the blue straggler region in the 
CMD. However, V107 is blended with a brighter star (see Figure[5j 
and therefore it is likely that the reference flux for this variable is 
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Figure 6. The period-luminosity relation in the V and / wavebands for the 
SX Phe stars in NGC 5024. The new SX Phe stars confirmed in this pa- 
per are plotted as solid red circles. The fil led and empty bl ack circles are 
the SX Phe stars discovered by All and ljeon et al.1 1 |2003|). respectively. 
Empty squares and triangles are the SX Phe stars from iDekanv & Kovfic j 
12009). The continuous straight line is the result of a least-squares fit to 
the stars pulsating in the fundamental mode. The short and long dashed 
lines correspond to the relations inferred for the first and second overtone 
pulsators, respectively, assuming the frequency ratios F/IO = 0.783 and 
F/20 = 0.571. The dotte d line in the top panel is the P-L relation calcu- 
lated bv ljeon et al] J2003I) from six SX Phe in NGC 5024. See Figure 12 in 
All for more details. 



over-estimated. Assuming that V107 is actually an SX Phe vari- 
able, which is an assumption consistent with its short period, then 
the blending could explain why it is brighter than the cluster P-L 
relation and why its colour places it outside of the blue straggler 
region. We therefore tentatively classify VI 07 as an SX Phe vari- 
able for which we cannot speculate on the mode(s) in which it is 
pulsating. 

In Sections [2.31 & [2~4l we tentatively confirm the variable na- 
ture of the candidates Wll and SX12 from SS11, but since our 
detection of the variability is border-line in these cases, we refrain 
from assigning these stars a V-number or speculating on their type 
of variability. Further follow-up observations of better precision 
will be required to properly confirm the variable nature of these 
particular stars. 
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Table 1. Properties of the new variable stars in NGC 5024. We list the mean V magnitude, mean V — I colour, detected frequency /, V-waveband photometric 
peak-to-peak amplitude A, and corresponding period P — f in columns 2, 3, 4, 5, and 6, respectively. The celestial coordinates at the epoch of our V reference 
image (heliocentric Julian date ~2455249.332 d) are listed in columns 7 and 8. The variable type is listed in column 9 along with the pulsational mode for the 
SX Phe variables in column 10. The numbers in parentheses indicate the uncertainty on the last decimal place. 
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V109 : P = 21.93 d 




Figure 7. Light curve from All and A12 for V109 phase-folded using a 
period of 21.93 d. V and / filter observations are plotted in black and red, 
respectively. The / filter light curve has been offset from the V light curve 
by 0. 1 mag for clarity. The fitted constant and sine-curve models are plotted 
as dashed and continuous curves, respectively. 



Finally in this Section, we report the discovery that the vari- 
able candidate VC10 from SS11 is actually a long period variable. 
Although we have ruled out the short term variations of ~0.3 mag 
amplitude claimed by SSI 1 (see their Figure 13), a careful inspec- 
tion of our corresponding light curve has revealed variations on the 
timescale of days with an amplitude of a few percent (discernable 
between observation nights in Figure [4). We have performed a pe- 
riod search in the range 10-1000 d again using PERI0D04 and we 
find the period with the strongest peak in the frequency spectrum is 
21.93±0.02 d. However, it is quite possible that the period we have 
found is an alias of the true period since our phase coverage at the 
timescale of days is quite low. 

In Figure [7] we plot the phase-folded light curve of VC10 at 
the period of 21.93 d using black and red points for the V and / 
photometric measurements, respectively. For clarity, we offset the 
/ light curve from the V light curve by 0.1 mag. Clearly the pho- 
tometric observations are well correlated between the two wave- 
bands. Furthermore, a sine-curve fit at the derived period of 21.93 d 
yields similar amplitudes of ~0.050 and 0.048 mag in the V and / 
wavebands, respectively, and Ax 2 values of ~ 150.6 and 169.7, re- 
spectively. Hence there is no doubt that this star is exhibiting low- 
amplitude long-period true variability. The star lies neatly on the 
red giant branch of the CMD, and we classify it as a semi-regular 
(SR) type variable. We assign the variable number V109 and report 
the variable star properties in Table [T] 

We comment that for V109, the variable nature of this star was 



not detected by SS11 for the simple reason that it is impossible to 
detect low-amplitude smooth variations at a 21.93 d period with 
only ~7.2 h of continuous observations. Therefore the discovery of 
VI 09 should be attributed to the work in this paper, as opposed to 
the discovery of V106-V108, which should be attributed to SS 1 1. 

In Table[2] available in full in electronic form (see Supporting 
Information), we provide our V and / time-series photometry for all 
the SX Phe and long-period variables in NGC 5024 from the works 
of A 1 1, A 12, and this papefj 



4 DISCUSSION 

SS11 consider a combination of several statistics to select their 74 
candidate variable stars. For ea ch of their ~10 4 light curves, they 
calculate an alarm statistic A jTamuz. Mazeh & North! l2006t) . an 
excess variance statistic Oxs iVaughan et al. 20031) , the Lomb pe - 
riodogram significance level F at the best-fit period jLomblll976h . 
and the RMS magnitude deviation r. They then devise the follow- 
ing combined criterion for selecting variable star candidates: 




(1) 



which selects 310 candidate variables. SS11 then rejected variable 
star candidates within 20 pix of the image edges and within 10 pix 
of known variable stars (although this last rejection criterion was 
obviously not applied properly), and used visual inspection to make 
the final selection of 74 candidates. 

As pointed out in A12, the alarm statistic A is flawed for 
variable star detection because of the normalisation by x 2 m its 
definition. Large amplitude variable stars have long runs of large 
residuals, but also a very high x 2 statistic, meaning that true large- 
amplitude variables yield undesirably small values for A. Therefore 
we are sceptical that the constraint on the A statistic applied in SS 1 1 
was actually useful for discriminating variable from non-variable 
stars. 

Both the Lomb periodogram significance level F and the ex- 
cess variance statistic Oxs are highly dependent on the accuracy of 
the estimated uncertainties on the photometric data points. Under- 
estimation of the uncertainties leads to smaller values for F and 
larger values for a X s, and vice versa. Since there is no discussion 



3 The light curves for all the RR Lyrae stars in NGC 5024 have already 
been published in Al 1 and A12. 
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Table 2. Time-series V and / photometry for all the SX Phe and long-period variables in NGC 5024 from the works of Al 1, A12, and this paper. The standard 
M st d and instrumental m ms magnitudes are listed in columns 4 and 5, respectively, corresponding to the variable star, filter, and epoch of mid-exposure listed in 
columns 1-3, respectively. The uncertainty on mj ns is listed in column 6, which also corresponds to the uncertainty on M st( j. For completeness, we also list the 
quantities f m f, f^if, and p from Equation 1 in All in columns 7, 9, and 11, along with the uncertainties o~ rc f and cyiff in columns 8 and 10. This is an extract 
from the full table, which is available with the electronic version of the article (see Supporting Information). 

Variable Filter HJD M std m im O m / ref f7 ref f m o~ di ff p 

Star ID (d) (mag) (mag) (mag) (ADUs^ 1 ) (ADUs"') (ADUs^ 1 ) (ADUs" 1 ) 

V67 V 2454940.17574 14.131 15.410 0.001 5984.296 1.689 1058.367 3.975 1.2186 

V67 V 2454940.19142 14.128 15.408 0.001 5984.296 1.689 1081.793 3.922 1.2230 

V67 / 2455219.49351 12.922 14.312 0.001 18675.923 3.782 131.660 12.956 0.7850 

V67 / 2455219.50345 12.919 14.310 0.003 18675.923 3.782 21.386 5.951 0.1031 



in SSI 1 of how the thresholds on F and Oxs were chosen, it is not 
clear how appropriate the adopted thresholds are. 

The RMS diagram presented in Figure 8 of SS11 shows 
that r > 0.01 mag for all stars fainter than ~ 17.5 mag in the R- 
waveband. Therefore the threshold on r chosen by SS11 does not 
discriminate variable stars from non-variable stars except for the 
brightest stars in the sample. When searching for variability via the 
r statistic, it is more appropriate to set a magnitude depedendent 
threshold determined from the RMS diagram itsel f (see the Vidrih 
index defined in Section 3.1 of lBramich et al.l2008t) since the RMS 
magnitude deviation is a magnitude dependent quantity that is gen- 
erally larger for fainter stars. 

Setting the variability detection threshold is a very important 
process that can be used to minimise the rate of false positive detec- 
tions and maximise the number of detections of true variables. Our 
investigation into the candidate variables presented by SS 1 1 clearly 
shows that the vast majority of their candidates are false positive de- 
tections with only a few detections of true variability, which means 
that they set their detection threshold too low. In fact, in SS 1 1, there 
is no assessment of the impact of random noise, the accuracy of the 
estimated photometric uncertainties, and the systematic trends in 
the light curves on the actual detection thresholds that should be 
applied in order to minimise the number of false positives. Assum- 
ing that the estimated photometric uncertainties are correct and that 
there are no systematic trends in the light curves, then Monte Carlo 
simulations employing the photometric uncertainties on the light 
curve data points may be used to determine the detection thresh- 
old required to achieve a desired rate of false positives. However, 
in practice, the estimated photometric uncertainties are usually in- 
accurate and systematic trends almost certainly exist in the light 
curves (as evidenced already in the light curve plots from both SS 1 1 
and this paper). Therefore, an appropriate detection threshold can 
only be determined from analysis of the distribution of the detec- 
tion statistic(s) for the light curves of the ensemble of stars under 
the assumption that the majority of stars are non- variable. 



5 CONCLUSIONS 

We have investigated the 74 candidate variable stars presented by 
SS11. We could not investigate four of these candidates because 
they are either too faint (three cases) or the DanDIA software failed 
to measure a reference flux on either reference image (one case). 
Of the remaining 70 candidates, we find that 14 are true variable 
stars, 10 of which were already discovered in previous works (al- 



though only 2 of these were available in the literature at the time of 
the SS11 paper submission), and four of which are new discover- 
ies. We assign the variable star numbers V106-V109 to the newly 
discovered variables. We classify V106-V108 as SX Phe type vari- 
ables whose discovery is attributed to SSI 1, and we classify V109 
as an SR type variable discovered in this paper. We also tentatively 
confirm the presence of true variability in two other candidates. 

From the analysis of our data, we strongly suggest that the 
other 54 candidate variable stars (or the majority of such) presented 
by SS 1 1 are false positive detections where systematic trends in the 
light curves have been mistaken for true variability. For 17 of these 
cases, we can explain the systematic trends as having been caused 
by their proximity to known RR Lyrae variables which are much 
brighter, and for 8 of these cases, the trends are due to their proxim- 
ity to saturated stars which reduce the local quality of the difference 
images. For the remaining 29 candidates, investigation of our light 
curves does not yield any evidence of true variability at the ampli- 
tudes reported by SS 1 1 , and any lower amplitude variations that we 
do see do not generally show the correlations between wavebands 
typical of true variability, which leads us to conclude that in these 
cases, the variations seen in the SSI 1 light curves are due to trends 
introduced during the reduction process. 
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APPENDIX A 

Briefly in this Appendix, we clarify some mistakes in the literature 
regarding variable star identification and astrometric coordinates. 

In Table 2 of All, there is a typo in the coordinates for 
V43. The coordinates should read: RA = 13 h 12 m 53.08 s , Dec. = 
+ 18°10'55.5", J2000. 

In All, we did not report the celestial coordinates of V53, 
which is highly blended. For completeness, we have measured 
these coordinates on the difference images where the variable star 
is not blended and we report them here as: RA = 13 h 12 m 55.85 s , 
Dec. = +18°10'35.5", J2000. These coordinates correspond to the 
epoch of our V reference image which is the heliocentric Julian 
date ~2455249.332 d. 
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We note that SS 1 1 arrive at the same conclusion as Al 1 about 
the previous misidentifications in the literature of the variables V57 
and V72, and the non-variability of V81, V82, and V83. 
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SUPPORTING INFORMATION 

Additional Supporting Information may be found in the online ver- 
sion of this article. 

Table 2. Time-series V and / photometry for all the SX Phe 
and long-period variables in NGC 5024 from the works of All, 
A 12, and this paper. 

Please note: Wiley-Blackwell are not responsible for the con- 
tent or functionality of any supporting materials supplied by the au- 
thors. Any queries (other than missing material) should be directed 
to the corresponding author for the article. 
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